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The correlation function of two identical pions interacting via Coulomb potential is 
computed for a general case of anisotropic particle's source of finite life time. The effect 
of halo is taken into account as an additional particle's source of large spatial extension. 
Due to the Coulomb interaction, the effect of halo is not limited to very small relative 
momenta but it influences the correlation function in a relatively large domain. The 
relativistic effects are discussed in detail and it is argued that the calculations have to be 
performed in the center-of-mass frame of particle's pair where the (nonrelativistic) wave 
function of particle's relative motion is meaningful. The Bowler-Sinyukov procedure to 
remove the Coulomb interaction is tested and it is shown to significantly underestimate 
the source's life time. 



1. Introduction 

The correlation functions of two particles with 'small' relative momenta provide 
information about space-time characteristics of parti cle's sources in high-energy 
nucleus-nucleus collisions, see the review ar ticlej™. Within the standard 'fem- 
toscopy' method, one obtains parameters of a particle's source, which is usually 
called the fireball, comparing the experimental correlation functions to the theoret- 
ical ones which are calculated in a given model. Such an analysis can be performed 
for pairs of non-identical or identical particles. In the former case, the correlation 
appears due to inter-particle interaction while in the latter one the interaction is 
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combined with the effects of quantum statistics. Since we usually deal with elec- 
trically charged particles, observed two-particle correlations are influenced, if not 
dominated, by the Coulomb interaction. The effect of the Coulomb force is treated 
as a correction and it is usually eliminate d fro m the experimental data by means 
of the so-called Bowler-Sinyukov procedur e! '^1^1 And then, the correlation function, 
which is obtained in such a way from experimental data, is compared with the 
theoretical correlation function of two non-interacting particles. The latter one is 
computed for a given source function - a distribution of the emission points - pa- 
rameterized usually in a gaussian form. The comparison with the experimental data 
provides parameters of the sourc e fun ction. 

Within the method of imagin^^l^, one obtains the source function not referring 
to its specific parametrization but directly inverting the functional dependence of 
the correlation function on the source function. The procedure of inversion involves 
the effect of quantum statistics as well as that of inter-particle interaction. The 
method provides essentially model independent information on the source space- 
time sizes but modeling is still needed to deduce the source life time which is 
coupled to the spatial size parameters. The method of one-dimensional imaging was 
successfully applied to experimental data- but the three-dimensional imaging is 
technically very complex and cumbersome . However, very first results of the three- 
dimensional imaging applied to NA49 data are presented in these proceedingJI^. 

The femtoscopy was applied to a large volume of experimental data on nucleus- 
nucleus collisions at SPS energjEl. The spatial size of particle's sources appeared 
to be comparable to the expected size a fireball created in nucleus-nucleus colli- 
sions while the emission time of particles was significantly shorter. It was predicted 
that at RHIC energies the emission time would be significantly longer due to the 
long lasting hydrodynar uic evo lution of the system created at the early stage of 
nucle us-nucleus c ollisionJ^^^. To a big surprise the experimental data obtained at 
show a very little, if any, change of the space-time characteristics 
of a fireball when compared to the SPS data. In particular, the emission time of 
particles appeared to be as short as 1 fm/c. This surprising result, which is now 
known as the 'HBT PuzzleH^HHI, 

demonstrates that cither we do not really un- 
derstand the dynamics of relativistic heavy-ion collisions or the principles of the 
femtoscopy are doubtful. In any case, the femtoscopy method need to be examined. 

In this presentation there is given a preliminary account of our study which is 
aimed to check three aspects of the standard method. Relativistic effects are, in 
our opinion, not satisfactory incorporated, as the non-rclativistic wave function of 
two interacting particles is often implicitly treated as a Lorentz scalar. It should be 
stressed that a transformation law of a wave function under Lorentz transformations 
is unknowrQ. However, the observed correlation functions are significantly different 
from unity for a small relative momenta when the relative motion of particles is 



^Using the Bethe-Salpeter equation, it has bee n recently shown that the hydrogen atom wave 
function experiences the Lorentz contractioJIS] under the Lorentz boost. 
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non-relativistic. Therefore, it is fully legitimate to use the non-relativistic wave 
function in the center-of-mass frame of two particles. However, it requires an explicit 
transformation of the source function to the center-of-mass frame. In Sec. S] the 
relativistic effects are discussed in detail and it is shown that a proper treatment 
of them leads to a numerically different result when compared with the standard 
procedure. 

The correlation function of two identical non-interacting bosons is expected to 
be equal to 2 for vanishing relative momentum of the two particles. The correlation 
function extracted from experimental data by means of the procedure, which is 
supposed to remove the Coulomb interaction, does not posses this property. The 
correlation function at zero relative momentum rather equals 1 + A with < A < 1. 
This fact can be explained referring to the concept of halcPQl. It assumes that 
only a fraction of observed particles, which equals A, comes from the fireball while 
the rest originates from the long living resonances. Then, we have two sources of 
particles: the fireball, which is rather small, and the halo with the radius given by 
the distance traveled by long living resonances. The complete correlation function, 
which includes particles from the fireball and the halo, equals 2 at exactly vanishing 
relative momentum. However, the quantum statistical correlation of two particles 
coming from the halo occurs at the relative momentum which is as small as the 
inverse radius of the halo. Since experimental momentum resolution is usually much 
poorer and such small relative momenta are not accessible, the correlation function 
is claimed to be less than 2 for effectively vanishing relative momentum. While the 
effect of halo is commonly believed to resolve the problem, the concept has been 
never examined for the particles which experience Coulomb interaction. In Sec. [6] 
we show that the effect of halo is not limited to very small relative momenta but due 
to the Coulomb repulsion the correlation function is modified for larger momenta. 

The Bowler-Sinyukov correction procedure, which is used to eliminate the 
Coulomb interaction from the experimental data, assumes that the Coulomb effects 
can be factorized out. The correction's factor is calculated for a particle's source 
which is spherically symmetric and has zero life time. We find such an approach 
rather inconsistent as the objective is to determine the source spatial parameters 
not assuming any source symmetry. The Bowler-Sinyukov procedure is examined 
in Sec. [7] wh ere we also comment on a similar test of the procedure performed in 
the papeJ^. 

We use the natural units, where c — h — 1, and the metric convention 

(+,-,-,-)■ 



2. Definition 

The correlation function C(pi,p2) of two particles with momenta pi and p2 is 
defined as 

dN 

(-(Pl,P2j - rfjv djV ' 
rfpi dp2 
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where ^pf^p^ and ^ is, respectively, the two- and one-particle momentum distri- 
bution. The correlation function can be written down in a Lorentz covariant form 



J- dpi ^ dp2 

where is the Lorentz invariant distribution. 

The covariant form ([T]) shows that the correlation function is Lorentz scalar 
which can be easily transformed from one to another reference frame. If the particle 
four-momenta, which are on mass-shall, transform as Pi — > p'j with i = 1, 2, the 
transformed correlation function equals 

C'(p'i,p'2) = C(pi(p'i),P2(p'2)). 

3. Nonrelativistic Koonin formula 

Within the Koonin modePSl, the correlation function C can be expressed in the 
source rest frame as 

C^(Pi,P2)= j d\^dhd?r2dt2D{viM)D{v2M)\^{T^'iy2)? , (2) 

where = -I- v^t;, '^{v'^, rj) is the wave function of the two particles and D{y, t) 
is the single-particle source function which gives the probability to emit the particle 
from the space-time point (i,r). The source function is normalized as 

j d^rdtD{r,t) = 1 . (3) 

After changing the variables r' ^ r, the correlation function can be written in 
the form 

C'(Pi,P2) = J dVidiidV2di2£'(ri -viti,ti)i?(r2-V2i2,<2)|*(ri,r2)|2 . 

Now, we introduce the center-of-mass coordinates 

r==r2-ri, R = p-(miri -f m2r2), 

t = t2-ti, T = ^{171111+112212), 

q = p-(m2Pi - miP2), P = pi+P2, 
where M = mi + 7122- Using the center-of-mass variables, one gets 

C(q)- J d^rdtDr{r~vt,t)\^^{r)\\ (4) 

where the 'relative' source function is defined as 

Dr{r,t)= /d3i?dTi^(R-^r,T-^t)i?(R+^r,r-|-^i). (5) 

J M M ' ^ M M ' ^ ' 
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To get Eq. (HI, the wave function was factorized as 

«'(ri,r2) = e'P^(^q(r) 

with </?q(r) being the wave function of the relative motion in the center-of-mass 
frame. Deriving Eq. ([4]), it has been assumed that the particles move with equal 
velocities i.e. vi = V2 = v which requires, strictly speaking, q = 0. However, one 
observes that |vi — V2I <C |vi| if |q| <C ^\pi\/mi where /x = mim^/M. Thus, the 
approximation vi = V2 holds for a sufficiently small center-of-mass momentum. 
We choose the gaussian form of the single-particle source function D{r,t) 



where r = {x, y, z) and the parameters r, Rx, Ry and Rz characterize the life time 
and sizes of the source. Specifically, the parameters r and R^ give, respectively, 



= (t^) EE J (frdtfD{v,t) , Rl = {x^) EE j d^r dtx^D{v,t) 



The relative source function computed from Eq. ([5]) with the single-particle 
source © is 



n I \ ^ x"^ t'^ 



IGn^RxRyRzT 



We note that the particle's masses, which are present in the definition ([5]), disappear 
completely in the formula This is the feature of the gaussian parameterization 
®- 

In the case of non-interacting identical bosons, the two-particle symmetrized 
wave function is 

*(ri,r2) = ^[e'Pi'-i+P^'-^ + e*P2ri-Hpir2j ^ -L[e'qr + g-iqrjgiPR^ 
V 2 V 2 

It gives the modulus square of the wave function of relative motion |iy9q(r)p = 
1 + cos (2qr) which in turn provides the correlation function equal to 

C(q) = 1 + exp [ - 4(T2(qv)2 + Rlql + R^^ + R^ql)] , (8) 

where q = (q^ , qy , qz) ■ We note that the 'cross terms' such as qxQz do not show up 
as the source function ([6]) obeys the mirror symmetry D{r,t) = D{—r,t). We also 
note that q often denotes the relative momentum pi — p2 not the center-of-mass 
momentum, which for equal mass nonrelativistic particles equals ^(pi — P2), and 
then, the factor 4 does not show up in the correlation function ([5]) of identical free 
bosons. However, we believe that using the center-of-mass momentum is physically 
better motivated. 
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4. Relativistic formulations 

There are two natural ways to 'relativize' the Koonin formula ([2]). The first one 
provides explicitly Lorentz covariant correlation function but it is applicable only 
for the non-interacting particles. The second one holds only in a specific reference 
frame but it is applicable for interacting particles as well. Below, we consider the 
two methods. We start, however, with the discussion of the Lorentz covariant form 
of the source function. 



4.1. Lorentz covariant source function 



The Lorentz covariant form of the gaussian parameterization of the source function 
(El) is 



D{x) 



VdetA 



1 



exp[ 



(9) 



47r2 """^ 2" 

where x^^ is the position four-vector and A^'' is the Lorentz tensor characterizing 
the source which in the source rest frame is 







1 

w, 












1 



(10) 



The source function as written in Eq. ([9]) obeys the normalization condition ([3]) not 
only for the diagonal matrix A but for non-diagonal as well. 

The source function ((9]) is evidently the Lorentz scalar that is 



VdctTV 



1 



VdetA 



^'(^') = ^^^^^Pi-n<^""'<] = ^4rFexp [-lx^A^^•'x,] = D{x) , 



47r^ 



^ ^m"^'' ^'^^ ^ L^^cr^'^^Lp" "^ith being the matrix of Lorentz 
transformation. We note that detA' = detL detA detL"^ — detA. 



where x',, 
msform 

The covariant relative source function ([7|) is given by 

Dr{x) 



VdetA , 
-1^ 2 '^^P 

lOTT-^ 



1 



(11) 



4.2. Explicitly covariant 'relativization' 

As follows from Eq. ([T]), the correlation function is a Lorentz scalar. Therefore, the 
Koonin formula Q can be 'relativized' demanding its Lorentz covariance. Let us 
write the formula 

C{pi,P2)^ j d^Xid\2D{xi)D{x2)\^{xi,X2)?, (12) 

where pi and Xi is, respectively, the four-momentum and four-position. Since the 
source function D{x) and the four-volume element dl^Xi are both the Lorentz scalars. 
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the whole formula ((T^) is covariant if the wave function ^{xi,X2) is covariant as 
well. In the case of non- interacting bosons the relativistic wave function '^(xi,X2) 
is 

^{xi,X2) = -l^le'P^^^+'P^''^ + e'P''''+'P''''). (13) 

As the function depends on the scalar products of two four-vectors, this is the 
Lorentz scalar. We note that the function ()13|) depends on two time arguments. 

Our further considerations are limited to pairs of identical particles and thus, 
we introduce the relative coordinates as 

X = X2 — Xi, X = + X2), 

(14) 

9=5(^1-^2), P^Pl+P2- 

We note that in the non-relativistic treatment the three-vectors r and q, which are 
given by the four- vectors x — {t, r) and q — {qo,<i), correspond to the inter-particle 
separation and the particle's momentum in the ccntcr-of-mass of the particle pair. 
This is, however, not the case in the relativistic domain. To get the center-of-mass 
variables, the four-vectors need to be Lorentz transformed. 
With the variables (fT4|) . the wave function (fT3|) equals 

*(x,X) = ^(e*«^ +e-*«^)e-*^^ , 
v2 

and the correlation function is found in the form 

C(g) = l + exp[-4q^(A^'')-igJ, 

which is explicitly Lorentz covariant. For the source matrix (|10p . the correlation 
function equals 

C{q) ^ 1 + exp [ - 4(<ZoV2 + qlRl + q^Rl + qlRl)] . (15) 

If |q| ^ |pi| with i = 1,2, then qo « qv, and the correlation function ([15]) exactly 
coincides with the non-relativistic expression ([8]) . This coincidence is not completely 
obvious as the time variables enter differently in the Koonin formula ([2]) and in the 
covariant one (H^). 

Let us consider the correlation function in the center-of-mass frame of the parti- 
cle's pair. We assume that the velocity of the center-of-mass frame in the source rest 
frame is along the axis x. Then, v = (v,0, 0) and qo « qv = q^v. The correlation 
function (jl5p . which holds in the source rest frame, equals 

C7(q) = 1 + exp [ - 4((^;V2 + Rl)ql + Rlql + Rlql)] . (16) 

As seen, the effective source radius in the direction x is y/R^^+v^r^ . We now 
transform the source function to the center-of-mass frame where the quantities are 
labeled with the index *. The center-of-mass source matrix (|10p . which is computed 
as 



Ar = L^A^'Lp^ 
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where 



7 —v"f 
-v"f 7 
10 
1 



with 7 = (1 — v^) equals 



Ar = 




















1 








1 



Then, the correlation function in the center-of-mass frame is found as 
= 1 + exp [ - A{j\v't' + Rl)ql^ + Rlqly + Rlql) 



(17) 



(18) 



As seen, the effective source radius along the direction of the velocity is elongated, 
not contracted as one can naively expect, by the factor 7. 



4.3. Non-covariant relativization 

The quantum mechanical description of two relativistic interacting particles faces 
serious difficulties. The problem is greatly simplified when the relative motion of two 
particles is non-relativistic (with the center-of-mass motion being fully relativistic). 
Since the correlation functions usually differ from unity only for small relative 
momenta of particles, it is reasonable to assume that the relative motion is non- 
relativistic. We further discuss the correlation functions taking into account the 
relativistic effects of motion of particles with respect to the source but the particle's 
relative motion is treated non-relativistically. In such a case, the wave function of 
relative motion is a solution of the non-relativistic Schrodinger equation. Thus, we 
compute the correlation function directly from the Koonin formula Q but we use it 
in the center-of-mass frame of the pair. For this reason we first transform the source 
function to this frame and then, after performing the integrations over xi and X2, 
we transform the whole correlation function, which is known to be a Lorentz scalar, 
to the source rest frame. 

As already stressed, we compute the correlation function in the center-of-mass 
frame of the pair and we use the relative variables (jl4p . Since the source function 
has the gaussian form ([9]) (with the non-diagonal matrix A), the integration over 
X can be easily performed and the correlation function equals 

C(q*) = / Dr{x,) |¥'q.(r^)|^ (19) 
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Fig. 1 . The Coulomb correlation function com- 
puted directly in the source rest frame ('fake 
covariance') and the correlation function com- 
puted in the center-of-mass frame of the pair 
and then transformed to the source rest frame 
('CM calculations'). 
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Fig. 2. The Coulomb correlation function 
C{qo, 0, 0) as a function of qo for the source pa- 
rameters: Rx = 2 fm, Ry = 2 fm, R^ = 3 fm, 
r = 2.5 fm and v = (0.9, 0, 0). 



where Dj.{x^) is the relative source function pip and Lp^^{v' ^) with r', = r» +v»t» 
is the non-relativistic wave function of relative motion. 

Although our aim is to compute the correlation functions of interacting particles, 
we start the discussion with the free identical bosons for the sake of comparison 
with the results of the previous section 14721 where the covariant 'relativisation' was 
presented. We again assume that v = (w, 0, 0), and then the source matrix is given 
by Eq. (fT7)l . The correlation function, which follows from Eq. ([TO]) , exactly coincides 
with the formula (|18p . To get the correlation function in the source rest frame, one 
performs the Lorentz transformation and obtains the formula p6p . Thus, the two 
ways of 'relativization' give the same result for non-interacting particles. This is 
not quite trivial as the time dependence of the Koonin formula ([TO)) and of the 
explicitly covariant one p2p is rather different. 

In the following sections we compute the correlation function of identical pions 
interacting due to Coulomb force, using the Koonin formula (fT9|) . Thus, we first 
compute the correlation function in the center-of-mass frame of the pair, and then, 
we transform it to the source rest frame. The correlation function of interacting 
particles is often calculated following the explicitly covariant method described in 
Sec. 14.21 And then, the non-relativistic wave function is treated as a Lorentz scalar 
function. Sometimes one argues that it must be a Lorentz scalar to guarantee that 
the right-hand side of Eq. \12\ is the Lorentz scalar. We believe that such a pro- 
cedure is incorrect and we refer to it as to 'fake covariance'. In principle, one can 
argue that the function '^{xi,X2) from Eq. (fT2| has to be a Lorentz scalar but it is 
unjustified to identify it with the non-relativistic wave function. The wave function 
is well defined in the center-of-mass of the pair and in the reference frames which 
move non-relativistically with respect to it. But properties of the wave function 
under Lorentz transformations are unknown, and thus, there is no reliable way to 
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Fig. 3. The Coulomb correlation function 
C(0, Qb, 0) as a function of qs for the source pa- 
rameters: Rx = 2 fm, Ry = 2 fm, Rz = 3 fm, 
T = 2.5 fm and v = (0.9, 0, 0). 



C(0,0,q|) 




q [MeV/c] 

Fig. 4. The Coulomb correlation function 
C(0, 0, gj) as a function of qi for the source pa- 
rameters: Rj: = 2 fm, Ry = 2 fm, ij^ = 3 fm, 
T = 2.5 fm and v = (0.9, 0, 0). 



transform it from one frame to another. We can also put it differently. The corre- 
lation function as defined by Eq. ([1]) is certainly a Lorentz scalar but there is no 
guarantee that the theoretical model ([2]) gives the correlation function which is a 
Lorentz scalar. However, we expect that the model works well in the center-of-mass 
frame of the pair where the non-relativistic wave function is well defined. Thus, we 
can compute the correlation function in this frame and then, we can transform it 
to an arbitrary frame, knowing that the correlation function is a Lorentz scalar. In 
this way, the non-covariant procedure circumvents the problem of unknown trans- 
formation properties of the non-relativistic wave function. 

Although the problem of transformation properties looks somewhat academic, it 
leads to a numerically significant effect. In Fig. [1] we show the correlation function 
of two identical charged pions computed in the center-of-mass frame of the pair 
and then transformed to the source rest frame. The function is compared to the 
correlation function which is directly computed in the source rest frame, treating 
the wave function as a Lorentz scalar. The source is assumed here to be spherically 
symmetric with i? = 2 fm and r = 0. The pair velocity with respect to the source 
chosen to be w = 0.9. (Details of the calculations of the Coulomb correlation func- 
tions are discussed in the next section.) As seen, the assumption of 'fake covariance' 
noticeably distorts the correlation function. 



5. Coulomb correlation functions 

In this section we compute, using Eq. (|19p . the correlation functions of two identical 
pions interaction via Columb potential. The calculations are performed for the 
anisotropic gaus siaii so urce of finite emission time ([9l fTO]) . We use the Bertsch- 
Pratt coordinateP^ES] ^^^^^ side, long. These are the Cartesian coordinates, where 
the direction long is chosen along the beam axis (z), the out is parallel to the 
component of the pair momentum P which is transverse to the beam. The last 
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Fig. 5. The Coulomb correlation function 
which includes the effect of halo for three val- 
ues of A equal 1.0 (upper line), 0.9 (middle 
line), 0.7 (lower line). 



'■'1 C{q) 
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Fig. 6. The free correlation function (upper 
line), the Coulomb correlation function (mid- 
dle line), the Bowler-Sinyukov correction fac- 
tor (lower line). 



direction - side - is along the vector product of the out and long versors. So, the 
vector q is decomposed into the Qq, Qs, and qi components. If the particle's velocity 
is chosen along the axis x, the out direction coincides with the direction x, the 
side direction with y and the long direction with z. We note that the correlation 
function of two identical free bosons in the Bcrtsch-Pratt coordinates in the source 
rest frame is 

C(q) = 1 + cxp [ - 4{qlRl + qlRl + qfRf) 



where Ro — + v'^t^, Rg = Ry and Ri ~ Rz- As seen, the source life time is 

mixed up with the size parameter R^- 

As well known, the Coulomb problem is exactly solvable within the non- 
relativistic quantum mechanics . The exact wave function of two non- identical 
particles interacting due to repulsive Coulomb force is given as 

^ (r) = W r(l + i^) f( ~ i^, 1, i{qr - qr)) , (20) 

where q = |q|, i] = fie^ /Sir with /i being the reduced mass of the two particles and 
±e is the charge of each of them; F denotes the hypergeometric confluent function. 
The wave function for the attractive interaction is obtained from (I20p by means 
of the substitution rj —rj. When one deals with identical particles, the wave 
function </?q(r) should be (anti-)symmetrized. The modulus of the symmetrized 
Coulomb wave function equals 



1 



V^{r)\' = -Giq) \F{~^^,l,^iqr-^v))\' + \F{-^-l,l,^{qr + qv))\' (21) 



q 



2Re(e*i'" ^(-«^, 1, Hqr - qr)) F*i-i^, 1, i{qr + qr))^ 
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Fig. 7. The correlation function C(go, 0,0) ex- 
tracted from the Coulomb correlation function 
by means of the Bowler-Sinyukov procedure 
(upper line) and the expected free correlation 
function (lower line). 



25 50 75 100 125 

q {MeV/c] 

Fig. 8. The correlation function C{0, qs^O) ex- 
tracted from the Coulomb correlation function 
by means of the Bowler-Sinyukov procedure 
(upper line) and the expected free correlation 
function (lower line). 



where G{q) is the so-called Gamov factor defined as 

= ^ ^ ■ (22) 

Substituting the modulus ([2T|) in Eq, (fT9|) . one finds the correlation function 
in the center-of-mass frame which is further transformed to the source rest frame. 
In Figs, m O and m we show the correlation functions C{qo, 0,0), C{0,qs,0) and 
C{0,0,qi), respectively, which are calculated for the following values of the source 
parameters Rx = 2 fm, Ry = 2 fm, R^ — 3 fm, t = 2,5 fm. The velocity of the 
particle's pair equals v = 0,9 and it is along the axis x. As seen, the correlation 
function C{qo, 0, 0) is qualitatively different than C(0, qs,0) and C(0, 0, qi). We note 
that except the Monte Carlo calculations presented in the papeJ^, where the source 
functions were generated according to the so-called blast-wave models, the results 
shown in Figs, [2l[3l [4] represent the first, as far as we know, relativistic calculations 
of the Coulomb correlation function for an anisotropic source of finite life-time. 

6. The Halo 

As mentioned in the introduction, the hald^Ql was introduced to explain the fact 
that, after removing the Coulomb effect, the experimentally measured correlation 
functions are smaller than 2 at vanishing relative momentum. The idea of halo 
assumes that only a fraction A of particles contributing to the correlation function 
comes from the fireball while the remaining fraction (1 — A) originates from long 
living resonances. Then, we have two sources of the particles: the small one - the 
fireball and the big one corresponding to the long living resonances. The source 
function has two contributions with A as a relative weight that is 

D{x) ^ XDf{x) + {l- X) Dh{x) , 
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Fig. 9. The correlation function C{0,0,qi) ex- 
tracted from the Coulomb correlation function 
by means of the Bowler-Sinyukov procedure 
(upper line) and the expected free correlation 
function (lower line). 
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Fig. 10. The correlation function extracted by 
means of the Bowler-Sinyukov procedure from 
the Coulomb correlation function, which in- 
cludes the halo, (upper line) and the expected 
free correlation function (lower line). 



where Df{x) and Dh{x) represent the fireball and halo, respectively. For non- 
interacting identical bosons, the correlation function is 

C(q) = 1 + + (1 - A)e-4-«'«' , (23) 

where both the fireball and halo are assumed to be spherically symmetric sources 
of zero life times; Rf and Rh are the radii of, respectively, the fireball and the 
halo. If Rh is so large that R~^^ is below an experimental resolution of the relative 
momentum g, the second term of the correlation function ([23| is effectively not 
seen, and one claims that C(q = 0) = 1 + A. 

We have included the halo in the calculations of the Coulomb correlation func- 
tions. The exemplary result is shown in Fig. [5] for three values of A: 1.0, 0.9, and 
0.7. For simplicity, the fireball and halo are spherically symmetric and have zero 
life times; Rf ^ 2 fm and Rh ~ 60 fm. 

7. The Bovirler-Sinyukov procedure 

As mentioned in the Introduction, the Coulomb effect is usually treated as a cor- 
rection and it is subtracted from the experimentally measured correlation functions 
by means of the Bowler-Sinyukov procedure. We first note that the Coulomb effect 
is far not small. In Fig. [6] we show the Coulomb and free correlation functions com- 
puted for the spherically symmetric source of zero life time with R = 2 fm. As seen, 
the correlation functions are qualitatively different from each other in the domain 
of small momenta q. Therefore, a method to subtract the Coulomb effect should be 
carefully tested. 

The Bowler-Sinyukov procedure assumes that the Coulomb effect can be factor- 
ized out, that is the correlation function can be expressed as 



C(q) = A{q) Cft.ec(q) , 



(24) 
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Fig. 11. The correlation function extracted by 
means of the 'diluted' Bowler-Sinyukov proce- 
dure from the Coulomb correlation function, 
which includes the halo, (upper line) and the 
expected free correlation function (lower line). 
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Fig. 12. The correlation function C(go, 0,0) 
extracted by means of the 'diluted' Bowler- 
Sinyukov procedure from the Coulomb corre- 
lation function, which includes the halo, (solid 
line) , the fitted free correlation function (dash- 
dotted line) and the expected free correlation 
function (dashed line). 



where Cfrcc(q) is the free correlation function and A{q) is the correction factor 
which depends only on q = |q|. 

The correction factor is actually the Coulomb correlation function which neglects 
the effect of quantum statistics i.e. the wave function is not symmetrized. The 
computation, which is, in particular, described in detail in the Appendix to the 
paper-"^, is performed in the ccntcr-of-mass frame where the particle's source is 
assumed to be symmetric and of zero live time. Thus, the correction factor is given 
by the formula 

A{q.)= J d\ Dr{r) \^g,{r)\^ ^ G{q,) J i?,(r) |^(-^, 1, i(g,r - q,r))|2, (25) 

where (pq_^ (r) is the Coulomb wave function pp]) and Dr (r) describes the spher- 
ically symmetric gaussian source of zero life time and of the 'effective' radius 
R = I \/i?o + i?s + i?i where i?o, Rs and Ri are the actual source radii; the Gamov 
factor G{q) is given by Eq. (|22p. Using the parabolic coordinates, the double inte- 
gration can be easily performed analytically in Eq. (|25p . and one is left with the 
one-dimensional integral which has to be taken numerically. Finally, the transfor- 
mation to the source rest frame should be performed. The exemplary result for 
i? = 2 fm is shown in Fig. [6l The pair velocity with respect to the source is as- 
sumed here to be so small that q = q^. As seen, the factor A{q) vanishes as q tends 
to zero. Thus, the (multiplicative) correction to Cficc(9) is infinite for q = 0. 

Once we are able to compute exact Coulomb correlation functions for an 
anisotropic source of finite life time, we can test whether the Bowler-Sinyukov 
procedure correctly subtracts the Coulomb effect for such a source and properly 
reproduces the free correlation functions. For this purpose, we first calculate the 
exact Coulomb correlation function in the pair center-of-mass frame, then we divide 



Towards reliable calculations of the correlation function 15 



C(0,q^,0) 

in- ^ 




25 50 75 100 125 

q [MeV/c] 

Fig. 13. The correlation function C{0,qs,0) 
extracted by means of the 'diluted' Bowler- 
Sinyukov procedure from the Coulomb corre- 
lation function, which includes the halo, (solid 
line) , the fitted free correlation function (dash- 
dotted line) and the expected free correlation 
function (dashed line). 
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Fig. 14. The correlation function C(0, 0, g;) 
extracted by means of the 'diluted' Bowler- 
Sinyukov procedure from the Coulomb corre- 
lation function, which includes the halo, (solid 
line) , the fitted free correlation function (dash- 
dotted line) and the expected free correlation 
function (dashed line). 



it by the correction factor A{q^,) computed according to Eq. (P5)) for a given effective 
source radius, and finally we transform the function to the source rest frame. To 
get source parameters, which are contained in the extracted 'free' correlation func- 
tion, the extracted function is fitted with the gaussian parameterization of the free 
function. The source parameters, which are used to compute the exact Coulomb 
correlation function, are called 'input parameters' while those which are obtained 
by fitting the extracted correlation function as the 'output parameters'. 

7.1. No halo 

We have first tested the Bowler-Sinyukov procedure for the case when all particles 
come from the fireball - there is no halo. Examples of the correlation functions 
extracted from the exact correlation functions by means of the Bowler-Sinyukov 



Table 1. The input and output source parame- 
ters, the radii and life time are given in fm. 





Ro 


Rs 


Ri 


T 


A 


input 


3.68 


2.50 


3.00 


3.00 


1.00 


output 


3.48 


2.49 


2.93 


2.69 


1.26 



Table 2. The input and output source parameters, the radii 
and life time are given in fm. 





Ro 


Rs 


Rl r 


A 


input 


4.39 


4.00 


6.00 2.00 


0.70 


output 


3.56 


3.60 


5.09 


0.51, 0.72, 0.64 



16 Maj & Mrowczynski 

procedure are shown in Figs. [3 [51 [HI The extracted functions are compared to the 
expected correlation functions of noninteracting bosons for the given source. The 
input and output parameters are given in Table [T] Actually, we have obtained three 
values of A as the extracted functions C{qo, 0, 0), C(0, Qs, 0) and C(0, 0, qi) are fitted 
independently from each other. However, all three output values of A are very close 
to each other, and we give only one value in the table. As seen, the correction 
procedure reproduces the source radii quite well. However, Ro is underestimated, 
and consequently, the source life time is reduced. We also note that the parameter 
A is increased by about 20%. 

7.2. Halo included 

The situation changes significantly when the halo is included in the consideration. 
In Fig. [10] we show the correlation function extracted from the Coulomb correlation 
function which includes the halo. In this illustrative example the fireball and halo 
are spherically symmetric and have zero life times; i?^^ = 2 fm and Rh = 60 fm. 
The correction factor is computed for R = 2 fm. For comparison we also show in 
Fig. [TOj the free correlation function computed for the double source of the fireball 
and halo with Rf = 2 im and Rh = 60 fm. We see that the extracted correlation 
function is badly distorted and the distortion extends far beyond R^^- 

The problem we face here is known and it can be partially resolved by extracting 
the 'free' correlation function not as in Eq. (f24|) but according to the formulaP^ 

C(q) = (l + X{A{q) - 1)) Cf,ec(q) • (26) 

The correction factor is 'diluted' that is only the contribution from the fireball is 
corrected. The correlation function extracted from the Coulomb correlation function 
by means of the formula (I26p is shown in Fig. (fTT|) together with the expected free 
function. The extracted function is less distorted than that one shown in Fig. [TOj 
but the distortion is still dramatic at small q and it extends beyond -R^^. 

We have applied the 'diluted' correction procedure to the Coulomb correlation 
function computed for an anisotropic source of finite life time with halo. Examples of 
the extracted correlation function are shown in Figs. [T^l [131 and [HI The extracted 
functions are presented together with the respective expected free functions and 
with the gaussian parameterizations of free function fitted to the extracted func- 
tions. When the latter functions were fitted with the gaussian parameterizations, 
the momenta smaller than Qmax, where ^max is the momentum for which the ex- 
tracted function has a maximum, were cut-off. The input and output parameters 



Table 3. The input and output source parameters, the radii 
and life time are given in fm. 





Ro 


Rs 


Rl r 


A 


input 


5.38 


4.00 


6.00 4.00 


0.70 


output 


4.07 


3.70 


5.08 1.88 


0.38, 0.60, 0.70 
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are given in Tabled The velocity of the particle's pair equals v = 0.9 and it is along 
the axis x. Since A is, as already mentioned, extracted separately from C((7o, 0,0), 
C(0, qs,0) and C(0, 0, qi), there are three output values of A which differ from each 
other. As also seen, all extracted source radii are significantly reduced. Since the 
output Ro is smaller than the output Rg, the output life time vanishes! A set of 
the input and output parameters with an extended input life time is presented in 
Table [3] The velocity of the particle's pair, as previously, equals v = 0.9 and it is 
along the axis x. We again observe a significant reduction of the source life time. 

An accuracy of the Bowler-Sinyukov procedure was also tested applying the 
Monte Carlo calculations presented in the papeJ^. And it was concluded that the 
source radii are reproduced quite well. However, the source functions were generated 
according to the so-called blast- wave models where the emission time is rather short. 
Then, as our findings also show, the Bowler-Sinyukov procedure works indeed very 
well. 

8. Conclusions 

We have presented here the preliminary account of our study of the two-particle 
correlation functions. We have argued that the calculations must be performed in 
the center-of-mass frame of the pair where a nonrelativistic wave function of the 
particle's relative motion is meaningful. We have computed the Coulomb correlation 
function of two pions coming from an anisotropic source of finite life time. The effect 
of halo has been also taken into account, and it has been shown that due to the 
Coulomb force the effect of halo extends for the particle's relative momenta far 
beyond the inverse halo radius. 

Having exact Coulomb correlation functions, the Bowler-Sinyukov procedure to 
remove Coulomb effect was tested. It was shown that the procedure works rather 
well when the halo is absent but with the halo the source radii are significantly 
reduced when compared to the original one. Since Ro is reduced more that Rg, the 
extracted life time of the source can be reduced even to zero. It might explain the 
'HBT puzzle' but a firm conclusion requires further systematic analysis which is 
still in progress. 
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